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We report evidence for the existence of a supercurrent of magnons in a magnon Bose-Einstein condensate pre¬ 
pared in a room temperature yttrium-iron-garnet magnetic film and subject to a thermal gradient. The magnon 
condensate is formed in a parametrically populated magnon gas, and its temporal evolution is studied by time-, 
frequency- and wavector-resolved Brillouin light scattering spectroscopy It has been found that local heating 
in the focal point of a probing laser beam enhances the temporal decrease in the density of the freely evolving 
magnon condensate after the termination of the pumping pulse, but it does not alter the relaxation dynamics of 
the gaseous magnon phase. This phenomenon is understood as the appearance of a magnon supercurrent within 
the condensate due to a temperature- and, consequently, magnetisation-gradient induced phase gradient in the 
condensate wave function. 


The term supercurrent is commonly assigned to a 
resistance-free charge current in a superconductor or to 
a viscosity-free particle current in superfluid 4 He or 3 He. 
More generally, this macroscopic quantum phenomenon 
relates to the collective motion of bosons (including real 
atoms [1, 2], Cooper pairs of fermions [3, 4] and bosonic 
quasiparticles [5-10]) coalescing into the same quantum state 
and, thus, constituting a specific state of matter described by 
a single, coherent wave function, a Bose-Einstein condensate 
(BEC) [11-14]. It is remarkable that any type of supercurrent 
is induced not by an external driving force but by a phase 
gradient in the BEC wave function. Up to now this prominent 
effect was only observed at rather low temperatures. How¬ 
ever, it is known that a BEC can be prepared even at room 
temperature by a strong increase in the density [15] of weakly 
interacting quasiparticles such as photons [16], plasmon- 
exciton polaritons [17] or magnons [18, 19]. Here, we report 
evidence for the existence of a supercurrent in a magnon 
BEC at room temperature in a ferrimagnetic insulator. The 
temporal evolution of the magnon BEC formed in a paramet¬ 
rically populated magnon gas was studied by Brillouin light 
scattering spectroscopy (BLS) [20, 21]. It has been found that 
heating in the focus of a probing laser beam locally enhances 
the decay of the magnon condensate but does not alter the 
relaxation dynamics of gaseous magnons. Our understanding 
of this effect is based on the temperature induced spatial 
variation in the saturation magnetization and, thus, on the 
variation in the local magnon frequencies across the heated 
sample area. Because the magnon condensate is coherent 
across the entire heated area, a spatially varying phase shift is 
imprinted into its wave function. This spatial phase gradient 
propels a magnon supercurrent that flows to the BEC areas 
around the focal spot and decreases, thus, the BEC density 
in the probing point. The detailed quantitative agreement 
between our experiments and a rate-equation model, which 
considers the temperature induced supercurrent jointly with 
a dispersion-caused magnon counter-current, provides a 
compelling proof for a supercurrent in a room-temperature 
Bose-Einstein magnon condensate. 


We have studied the appearance of a magnon supercurrent 
by observing the temporal evolution of a magnon BEC in an 
in-plane magnetized single crystal yttrium iron garnet (YIG, 
YsFe 50 i 2 ) film. YIG is a ferrimagnetic dielectric material 
with unique low magnon damping [22] allowing both for ef¬ 
ficient thermalization of externally injected magnons and for 
a several hundred nanosecond-long BEC lifetime. In order to 
achieve Bose-Einstein condensation [19, 23], magnons are in- 
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FIG. 1. Magnon spectrum of a 5.6 /xm thick YIG film magnetized in¬ 
plane by a bias magnetic field H = 1690 Oe shown for the wavevec- 
tor q perpendicular (left part) and parallel (right part) to the applied 
field. For both wavevector directions the first 31 thickness modes are 
shown. The red arrow illustrates the magnon injection process by 
means of parallel parametric pumping. N p - total number of para¬ 
metrically excited magnons at cj p /2; N c - number of BEC magnons 
at lj c — oxmin; Nb - number of gaseous magnons near u; m in and q m i n ; 
N g - number of magnons in the parametrically overpopulated gas of 
magnons below uj p / 2. The inset schematically shows the geometry of 
the in-plane magnon wavevector q, the bias field H, and the pumping 
microwave magnetic field. 
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FIG. 2. Experimental setup, a, Schematic illustration of the experimental setup. In the upper part of the figure the microwave circuit consisting 
of a microwave source, a switch and an amplifier is shown. This circuit drives a microstrip resonator, which is placed below the in-plane 
magnetized YIG film. Light from a solid-state laser (A = 532 nm) is chopped by an acousto-optic modulator (AOM) and guided to the YIG 
film. There it is scattered inelastically from magnons, and the frequency shifted component of the scattered light is selected by the tandem 
Fabry-Perot interferometer, detected, and analyzed in time. The time diagram b presents the relative timing of the microwave pump pulse, the 
laser pulse, and the detected BLS signal, c, Infrared photo of the YIG sample showing the hot spot in the focus of the probing laser beam, d, 
Schematic illustration of a supercurrent flowing in the magnon BEC subject to a thermal gradient. Red arrows - outflow of the thermally induced 
supercurrent J T from the hot focal spot. Blue arrows - contraflow of the dispersive current Jdi s - w c and cc® - magnon BEC phases in the hot 
and in the cold parts of the sample, respectively. Sip- temperature induced phase accumulation in the magnon BEC; R - radius of the hot spot. 


jected into the YIG spin system via parallel parametric pump¬ 
ing [24], which currently is considered to be the most efficient 
technique for magnon excitation over a large wavevector range. 
The process can be described by a splitting of a photon of a 
pumping electromagnetic wave with nearly zero wavevector 
and the pumping frequency c into two magnons with op¬ 
posite wavevectors ±q and frequency cj p / 2. This magnon in¬ 
jection process is schematically illustrated by the red arrow in 
Fig. 1, which presents the magnon spectrum calculated for the 
given experimental conditions. 

The strength of the bias magnetic field H = 1690 Oe is cho¬ 
sen to allow for magnon injection slightly above the ferromag¬ 
netic resonance frequency cu FMR , where the parallel pump¬ 
ing process has its highest efficiency [25]. When the injected 
magnons thermalize through scattering processes conserving 
both their number and the total energy [26, 27], a BEC may 
be formed at the lowest energy state of the energy-momentum 
spectrum even at room temperature [19, 23]. 

The experimental setup, which consists of a YIG film, a 
microwave circuit and a BLS system, is schematically shown 
in Fig. 2a. The YIG sample is 5 mm long and 1mm wide. 
The single-crystal YIG film of 5.6 fim thickness is grown in 
the (111) crystallographic plane on a gadolinium gallium gar¬ 
net (GdsGa 50 i 2 ) substrate by liquid-phase epitaxy. The mi¬ 
crowave pumping circuit comprises a microwave source, a 
switch, and an amplifier. This circuit drives a 50 /tm wide mi¬ 
crostrip resonator, which is placed below the YIG film. The 
resonator concentrates the applied microwave energy and in¬ 
duces a microwave Oersted field oriented along the bias mag¬ 
netic field H , thus realizing conditions for the parallel para¬ 
metric pumping mechanism. The pumping pulses of r p = 2 /is 
duration are applied with a repetition time of 1 ms allowing for 
magnetic and temperature equilibration of the YIG film. 

In our experiment, a focused laser beam combines the 
role of the magnon probe in the BLS experiment with the 


role of the local sample heater. The heating time is adjusted 
by chopping the probing laser beam using an acousto-optic 
modulator (AOM). The laser pulse duration r L is varied 
between 6 /is and 80 /is. By changing r L the setup allows one 
to change the heating time interval before the application of 
the microwave pumping pulse and, consequently, to control 
the sample temperature in the probing point (see Fig. 2b). The 
modulated probing beam is focused onto the surface of the 
YIG film sample in the middle of the microstrip resonator 
(see Fig. 2a and 2c), where it has a maximum peak power of 
9.5 mW. The focal spot has a radius R « 10 (im and, thus, is 
approximately twice smaller than the 50 fi m wide parametric 
pumping area. The scattered light is deflected by the beam 
splitter to the tandem Fabry-Perot interferometer, which 
selects the frequency shifted component of the scattered 
light. The intensity of the inelastically scattered light, which 
is proportional to the magnon density in the probing point, is 
detected and analyzed in time. 

In order to understand the of the magnon BEC, one first 
needs to separate the effects caused by the spatially uniform 
change of the sample temperature from those caused by the 
formation of a temperature gradient. Therefore, in a first 
experiment, we combined BLS probing at low laser power 
with an uniform heating of the YIG sample by a hot air flow. 
The sample temperature was measured by an infrared camera. 
Magnons were detected near their spectral energy minimum 
in a frequency band of 100 MHz spectral width. 

Figure 3a shows the typical dynamics of the magnon BEC 
in this case. During the action of the pump pulse the magnon 
density increases and saturates. After the parametric pumping 
is switched off, the magnon density jumps up due to the 
previously reported intensification of the BEC formation in 
the freely evolving magnon gas [23]. Afterwards the magnon 
density exponentially decreases due to the conventional 
spin-lattice relaxation mechanisms. This behavior is common 
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for all temperatures in our experiment. Some decrease in the 
steady-state magnon density observed at higher temperatures 
can be related both to higher magnon damping [22] and to 
lower efficiency of the parametric pumping in this case. The 
latter effect is straightforward: Due to a lowered saturation 
magnetization at higher temperatures the parametrically 
populated spectral magnon branch shifts towards lower 
frequencies, see Fig. 3b. This shift increases the wavenumber 
of the injected magnons and, consequently, the threshold 
power of the parametric generation [25]. 

Now we focus on the temperature-gradient dependent be¬ 
havior of the magnon condensate. Figure 4 shows the evo¬ 
lution of the magnon density at the bottom of the spin-wave 
spectrum for four different heating times, i.e. for four values 
of the temperature gradient, at a laser power of P L =9.5 mW. 
For comparison, the black curve presents the time-dependent 
magnon dynamics measured using a much lower laser power 
of 0.4 mW. Similar to Fig. 3a, the BLS signal rises sharply after 
the microwave pumping pulse is switched off due to the in¬ 
tensification of the BEC formation process [23]. Afterwards, 
the magnon density decreases. In the case of low-power laser 
probing this decrease has an exponential form with a charac¬ 
teristic decay time r ~ 240 ns, which corresponds well to the 
conventional values of a linear magnetic decay in YIG films 
[22] in the vicinity of cj m i n > i.e. for both the magnon BEC and 
the gaseous magnons near the bottom of the spectrum. 

At higher laser power the observed decrease in the magnon 
density cannot be described by one single exponential func¬ 
tion anymore. Two different regimes of the magnon decay are 
clearly visible from the slopes of the BLS signals presented in 
Fig. 4. In the first regime, at high BLS intensities, the decay of 
the magnon density depends on the heating duration, i.e. on 
the temperature gradient in the laser focal spot. In the second 
regime, at low BLS intensities, this decay is temperature 
independent and its rate coincides with the one observed 
in the low-laser-power experiment. It is remarkable that for 
all heating durations the transition between these two decay 
regimes occurs at approximately the same magnon density 
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FIG. 3. Temperature-dependent temporal dynamics of the magnon 
BEC. a, Time-resolved BLS intensity for different uniform tempera¬ 
tures of the YIG film, b, Influence of the YIG film temperature on the 
spin-wave dispersion (same color code for temperature as in a. 



FIG. 4. Temporal dynamics of the magnon BEC under local laser 
heating. The increase of the laser pulse duration, and thus the heat¬ 
ing time, leads to a faster decay of the magnon BEC in the measure¬ 
ment spot, i.e. laser spot. The relaxation time of the residual gaseous 
magnons at the bottom of the spin-wave spectrum is independent of 
the heating. The parametrically injected magnons at cu p /2 show no 
temperature related change in their nonlinear relaxation dynamics af¬ 
ter the end of the pumping pulse. 


level, which can be associated with the transition from the 
condensed to the gaseous magnon phase. 

It is also worth to analyze the evolution of parametrically 
injected magnons (see Fig. 4), which play the important role 
of a magnon source for all other spectral states of the magnon 
gas and, thus, can strongly influence the BEC dynamics. In 
the beginning, this dense magnon group rapidly decays due 
to intense four-magnon scattering processes, see e.g. Ref. [28]. 
This decay closely relates to the jump in the BEC population 
number shown in Fig. 4. In the course of time, the density of 
the parametrically excited magnons decreases, four-magnon 
scattering decelerates and, eventually, the magnon relaxation 
attains to its intrinsic value corresponding to the decay rate of 
the gaseous magnons at the bottom of the magnon spectrum. 
As is clearly seen in Fig. 4, this nonlinear relaxation behaviour 
is independent of the laser heating and, thus, cannot be the 
reason for the temperature-gradient-related BEC dynamics. 

We now present a model to describe our experimental find¬ 
ings. The magnon BEC is a spontaneously established co¬ 
herent ground state, which possesses well defined frequency, 
wavevector, and phase [29]. The local laser heating process 
locally changes the saturation magnetization and, thus, in¬ 
duces a weak frequency shift Suj c between different parts of the 
magnon condensate. In the course of time this frequency shift 
results in an increasing phase gradient in the magnon BEC. As 
a result, a phase-gradient-induced magnon current or, in other 
words, a magnon supercurrent, flowing out of the hot region 
of the focal spot, is excited. This efflux reduces the density of 
the magnon BEC in the probing point. After some time, the 
enhanced decrease in the magnon density results in the disap- 
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pearance of the condensate, and thus in the disappearance of 
the supercurrent. Consequently, this leads to the recovery of 
the conventional relaxation dynamics associated with a resid¬ 
ual incoherent gaseous magnon phase. 

We would like to emphasise that the laser heating process 
decreases locally the saturation magnetization via thermal ex¬ 
citation of high-energy magnons with terahertz frequencies 
but practically does not increase the local population of low- 
energy magnons. As a result, the number of thermal magnons 
in the low-energy spectral region remains negligibly small in 
comparison to the number of magnons originating from the 
parametric pumping process and cannot visibly affect the BEC 
dynamics. Moreover, the observed temperature-independent 
relaxation dynamics of the incoherent gaseous magnon phase 
(see Fig. 4) provides unambiguous evidence for the fact, that 
the interaction between high- and low-energy magnons is neg¬ 
ligibly small: within experimental resolution a diffusively ex¬ 
panding cloud of thermally excited high-energy magnons does 
not drag the low-energy gaseous magnons out of the heated 
spot and, therefore, does not contribute to the enhanced local 
decay of the magnon BEC. 

To quantify our interpretation, that a magnon supercurrent 
is acting, recall that a BEC of magnons can be described by 
the Gross-Pitaevskii equation for the complex wave function 
C(r,t) = V^c ex P(^)> where N c (r,t) and p{r, t) are the 
magnon number density per unit volume and the BEC phase, 
respectively. This equation is the well known equation for the 
envelope of a narrow wave packet (see e.g. Refs. [28, 30]). 
It conserves the total number of condensed magnons Af c = 
f N c (r , t) dr and leads to the continuity equation for the BEC 
density 7V c (r, t): 


dN c 

dt 


+ V-/ = 0, 


N c ^ d 2 uj(q) dp 
2 “ dqidqj drj 


(1) 


Here /(r, t) is a supercurrent formed by the BEC magnons. In 
the (x, y) -plane of a YIG film (see Fig. 2d) it has two compo¬ 
nents 


J x = N c D x dip/dx, J y = N c D v d(p/dy , (2) 

where D x = d 2 ui(q)/(2 dq x ) and D v = d 2 ui(q )/(2 dq 2 ) are 

the anisotropic dispersion coefficients. At our experimental 
conditions D x ~ 21 D y and, thus, J x J y . This allows us to 
neglect J y in the global balance of the magnon numbers and to 
simplify the problem to a ID case along the x-axis, see Fig. 2d. 

There are two reasons for the x -dependence of the BEC 
phase p in our experiment. The first one is the already men¬ 
tioned temperature dependence of uj c . Within the hot spot of 
radius R centered at x = 0 (i.e. for \x\ < R) the tempera¬ 
ture T{pc) is higher than the temperature T 0 of the rest of the 
film, (see Fig. 2d). Since in an in-plane magnetized YIG film 
duo c ( T)/dT <0, the BEC frequency in the spot is smaller than 
outside: 8uj c {x) = u c (T(x)) — cj c (Xb) < 0. Correspondingly, 
the phase accumulation Sp(x) = 5w c (x)t inside of the spot 
is smaller than in the surrounding cold film. Therefore, the 
phase gradient dSp(x)/dx is positive for x > 0 and negative 


for x < 0. It means that a thermally induced supercurrent 
flows out from the spot (mostly in x -direction) as it is shown 
by red arrows in Fig. 2d: 

J T = sign (x)N c D x • (3) 

This outflow decreases the magnon BEC density N c (x ) in the 
spot, \x\ < R , with respect of that in the cold film, where 

N c {x > R) = 7V C °. 

Spatial deviations in the density N c (x) of the magnon con¬ 
densate constitute the second reason for the variation of its 
phase dp/dx ^ 0. It results in an additional contribution 
to the supercurrent given by Eq. (1), which can be named a 
“dispersive” current Jdi s . In order to estimate J^[ s notice that 
the ID Gross-Pitaevskii equation has the self-similar solution 
C(x,t) = J-'(^) with £ = x 2 /(D x t), where the function 
R(£) satisfies the ordinary differential equation t^dR/dt; = 
i(dF/d £ + 2 ^d 2 R/d^ 2 ). It describes the well known phe¬ 
nomenon of dispersive spreading of a wave packet with the 
width 5x = y/D x t permanently increasing in time. It is worth 
to notice that exactly the same law is satisfied by the diffusion 
process with the diffusion current Jdif = D x (dN c {x)/dx). 
Ignoring the difference between the self-similar profiles in the 
dispersion and the diffusion processes (it is well below the res¬ 
olution of our experiment), we can estimate the dispersive cur¬ 
rent Jdis by its diffusion counterpart: 

Jdis 0*0 - -sign(ar) D x . (4) 

This current is directed to the lower N(x) region, i.e. towards 
the hot spot as shown in Fig. 2d by blue arrows. 

A consistent description of the evolution of a parametrically 
driven magnon system towards BEC may be achieved in the 
framework of the weak wave turbulence theory [31], which 
was further developed in Ref. [28] to describe the spin-wave 
turbulence under parametric excitation. The main tool of this 
theory is a kinetic equation for the occupation numbers n(q) 
of waves (magnons in our case). However, this ambitious aim 
can only be solved in the future. Here we restrict ourselves 
to the analysis of a crucially simplified version of the kinetic 
equations - by rate equations for the total number of magnons 
in particularly chosen parts of the g-space shown in Fig. 1. The 
first part is the BEC with density N c (orange dot in Fig. 1). The 
second one is the magnon gas with density N b at the bottom 
of the magnon spectrum in close vicinity to the BEC (blue area 
Mb surrounding the orange BEC point in Fig. 1). These parts 
are directly coupled via four-magnon scattering processes 
(Mb BEC) and are simultaneously detected due to the 
finite resolution of our frequency- and wavenumber-resolved 
BLS setup. The third part is given by the magnon gas area M g 
above Mb. In our simplified model the parametrically injected 
magnons of density N p (magenta dot in Fig 1) populate first 
the gas area M g (7V P —» M g ) and afterwards move to the 
bottom part Mb of the magnon spectrum due to four-magnon 
scattering processes between these areas (M g ^ Mb). 
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In the rate equations for the magnon numbers N g , N b and 
N c we have to account that the kinetic equations for the four- 
magnon scattering processes conserve the total number of 
magnons and that the leading term for the magnon flux from 
the j to the i sub-system is proportional to Nj and, thus, can 
be written as AijNJ with Aij being dimensional (s _1 ) phe¬ 
nomenological constants. 

Thus, the rate equations take the following form: 

IT = - r « + - A*" 3 ' + ’ < 5 *> 

dN h , , 

-fir- = -r b N b + A gh Nl - A he N b 3 

-A bc (N* - N* r )G(N h - N CI ) , (5b) 

^ = -r c N c + A bc (N g - iV c 3 r )0(AT b - N CI ) - (5c) 

Here T g , I\> and T c are relaxation frequencies of correspond¬ 
ing magnons in the processes that do not conserve magnon 
numbers (mainly caused by spin-orbit interaction). Their ex¬ 
perimental values are close enough and for simplicity are re¬ 
placed by the same mean value T = 6.1 MHz. 

The second term in the right-hand-side of Eq. (5a) repre¬ 
sents the external magnon source and is proportional to the 
number of parametrically pumped magnons N p . The factor 
exp(—T o t) models the processes of magnon thermalization af¬ 
ter the end of the pump pulse. The terms oc A^ in Eqs. (5) 
describe the magnon fluxes between the chosen spectral ar¬ 
eas, leading finally to the population of the BEC state. The 
flux Mb M c contains the Heaviside function ©(TVb — N cr ) 
that involves 7V cr —a critical number of magnons at which the 
chemical potential /i of the magnon gas reaches The 

function ©(iVb — N cr ) = 1 for — N cr > 0 and 0 oth¬ 
erwise. Therefore, for Nb < N cr there is no flux of magnons 
to the condensate. For Nb > N cr this flux appears and an ex¬ 
cess of magnons Nb — N cr populates the condensate. 

As we have discussed above, the supercurrent J in Eq. (5c) 
consists of two parts: J = J T + J^ s . Their modelling is, 
probably, the most delicate issue. We adopt the simplest 
version, assuming a constant temperature gradient from the 
center (x = 0) of the heated spot to R : 

J = —B s N c 6w c t + B dis (N° - N c ). (6) 

Here, for the coefficients we have B s ~ is ^ —D x /R. 

The results of the numerical solution of Eqs. (5) and (6) for 
different T in the center of the hot spot and, accordingly, for 
different 6w c (T), are shown in Fig. 5 by color solid lines to¬ 
gether with the respective experimental data. One clearly sees 
the good agreement between the theoretical and experimental 
curves. The initial rise of the theoretical curves occurs faster 
than in the experiment. This is likely due to our simplified 
model with only three zones in the g-space. We ignore here 
the multistage character of the N p M g magnon transfer. 

In a numerical solution of our model we have choosen the 
initial BEC value 7V C = 0 because the gaseous magnons do 
not populate the lowest energy state in the case of powerful 



FIG. 5. Theoretically calculated magnon dynamics in a thermal 
gradient. Theoretical dependencies (colour lines) of the observable 
magnon densities ( N c + Nb) calculated according to Eqs. (5) for dif¬ 
ferent temperature dependent shifts 6uj c {T) of the BEC frequency uj c 
are shown in comparison with the corresponding experimental data 
(gray curves). The inset shows the calculated temporal behaviour of 
iVc, AT, and N g for the case of the strongest heating of the magnon 
BEC. 

pumping [23]. In such cases, no magnon condensation occurs 
during the pump action and the BEC in the global energy mini¬ 
mum (lj m in, &min) is formed after the termination of the pump 
pulse. The initial value of Nb = 135 (in units of BLS counts) 
was taken from our experiment (measured value TVb + N c at 
t = 0 in Fig. 5). The initial values N p = 250 and N g = 1000 
are found from the best fit to the reference case of no heating, 
shown by a black line in Fig. 5. 

In the comparison of our model with experiments per¬ 
formed at different spot temperatures T, we used only one 
T-dependent parameter Slj c (T ), related to the temperature 
dependence of the saturation magnetization M(T ) of YIG 
[22]: at standard ambient conditions dM(T)/dT ~ —4G/K 
giving duj c /dT = — 116 kHz/K. For the largest frequency shift 
(Scu c = 550 kHz), assumed in our calculations, it requires 
a rather mild heating of the spot by about 4.7 K. This value 
agrees well with a numerically found temperature rise of 5.7 K 
in the center of the focal spot. This value was obtained by solv¬ 
ing a 3D heat conduction model in the COMSOL Multiphysics 
simulation package using known values of the light absorption 
coefficient, heat capacity, duration of the laser pulse, etc. 

To summarize, we used Brillouin light scattering spec¬ 
troscopy to study the time-dependent evolution of the magnon 
density in the narrow frequency and wavevector vicinity of the 
bottom of the magnon spectrum. We have found that local 
laser heating of the probing spot visibly enhances the tempo¬ 
ral decrease in the population of the magnon Bose-Einstein 
condensate freely evolving after the termination of paramet¬ 
ric pumping. At the same time, the relaxation dynamics of 
both gaseous and parametrically injected magnons was not 
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changed. In order to interpret the experimental results we pro¬ 
posed and analyzed a model based on rate equations for the 
gaseous and condensed magnons, taking into account both the 
temperature induced magnon supercurrent and the diffusion- 
like dispersive magnon counterflow via the boundary of the 
hot spot. In these equations, the position of the magnon en¬ 
ergy minima hw m i n (T) was assumed as the only tempera¬ 
ture dependent parameter. The surprisingly high quantita¬ 
tive agreement between the measured data and the analyti¬ 
cally modeled magnon evolution at the bottom of the magnon 
spectrum serves as strong support in favor of the existence 
of a temperature-induced supercurrent in the Bose-Einstein 
magnon condensate. Moreover, the occurrence of the super¬ 
current directly confirms the phase coherency of the observed 
magnon condensate and fuels expectations of the realization of 
phase-controlled magnon transport [32,33], which can be used 
in future magnonic devices for low-loss information transfer 
and processing [34]. Therefore, with this demonstration of a 
magnon supercurrent in a magnon BEC we have opened the 
door to studies in the general field of magnonic macroscopic 
quantum transport phenomena at room temperature as a novel 
approach in the field of information processing technology. 
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